Highly Anisotropic Gap Function in Borocarbide Superconductor LuNi2B2C 
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The thermal conductivity of borocarbide superconductor LuNi2B2C was measured down to 70 mK 
(Tc/200) in a magnetic field perpendicular to the heat current from iif = to above IIc2 = 7 T. As 
soon as vortices enter the sample, the conduction at T grows rapidly, showing unambiguously 
that delocalized quasiparticles are present at the lowest energies. The field dependence is very 
similar to that of UPta, a heavy-fermion superconductor with a line of nodes in the gap, and very 
different from the exponential dependence characteristic of s-wave superconductors. This is strong 
evidence for a highly anisotropic gap function in LuNi2B2C, possibly with nodes. 

PACS numbers; 74.70. Dd, 74.25.Fy, 74.60.Ec 



The vast majority of known superconductors are char- 
acterized by an order parameter with s-wave symmetry 
and a gap function which is largely isotropic and without 
nodes (zeros). Only four families of materials are seri- 
ously thought to exhibit a superconducting state with 
a different symmetry: (1) heavy-fermion materials, such 
as UPta where a line of nodes in the gap function has 
clearly been identified (2) high- Tc cuprates, such 
as YBa2Cu307 where the order parameter was clearly 
shown to have d-wave symmetry [H; (3) the ruthenate 
Sr2Ru04, where there is strong evidence for a triplet or- 
der parameter and (4) organic conductors, such as k- 
(ET)2Cu[N(CN)2]Cl where there is growing evidence for 
unconventional superconductivity A major outstand- 
ing question is the nature of the microscopic mechanism 
responsible for superconductivity in any of these materi- 
als. The unconventional symmetry of the order parame- 
ter is evidence for a pairing caused by purely electronic 
interactions and not mediated by phonons. For example, 
the proximity to magnetic order which is found in all four 
families of superconductors has led to the suggestion that 
spin fluctuations are responsible for Cooper pairing, as is 
thought to be the case in superfluid ^He. 

The presence of nodes in the gap function is gener- 
ally associated with unconventional (non-s-wave) sym- 
metries. These nodes are typically inferred from the ob- 
servation of quasiparticle excitations at energies much 
lower than the gap maximum Aq, as reflected for exam- 
ple in the power law temperature dependence of various 
physical properties, such as London penetration depth 
and ultrasonic attenuation at T <C Tc. Another way of 
detecting low-energy quasiparticles is to excite them by 
applying a magnetic field which introduces vortices in the 
material, so that the superfluid flow around each vortex 
Doppler shifts the quasiparticle energy. In certain limits, 
the quasiparticle response is the same whether induced by 
a thermal energy ksT or by a field energy ~ Aq ^/BJB^, 



where Bc2 — Hc2, the upper critical field ||]. 

In this Letter, we turn our attention to another class 
of superconductors: the borocarbides LNi2B2C (where L 
= Y, Lu, Tm, Er, Ho, and Dy) Q. It has generally been 
thought that these materials are described by an order 
parameter with s-wave symmetry and pairing which pro- 
ceeds via the electron-phonon coupling [0-^. However, 
there is recent evidence for low-energy excitations in the 
superconducting state, whether from the anomalous field 
dependence of the specific heat [|l^,|ll|l and the microwave 
surface impedance [ pdp^ ], or from the presence of scat- 
tering below the gap in Raman measurements [|l3| . This 
has been interpreted in terms of an anisotropic s-wave 
gap pO 11| (see also reference H). 



Here we present compelling evidence that the gap func- 
tion of LuNi2B2C is highly anisotropic, with a gap min- 
imum Aniin at least 10 times smaller than the gap max- 
imum, A„iin < Aq/IO, and possibly going to zero at 
nodes. This statement is based on the observation of 
delocalized quasiparticles at very low energies, as mea- 
sured directly by heat transport. Indeed, quasiparti- 
cle conduction is induced by a magnetic field as low as 
Hci — -ffc2/100 and it grows linearly with field, in dra- 
matic contrast with the exponentially activated transport 
seen in Nb, for example, where it results from tunneling 
between the localized states bound to the core of adja- 
cent vortices. Such pronounced anisotropy challenges the 
current view on the nature of superconductivity in boro- 
carbides. It suggests either a new family of unconven- 
tional superconductors (with symmetry-imposed nodes 
in the gap function) or s-wave superconductors with more 
anisotropy than has ever been seen before. In either case, 
the role of phonons as a proposed pairing mechanism may 
have to be critically re-examined. 

LuNi2B2C is an extreme type-II superconductor and a 
non-magnetic member of the borocarbide family with a 
superconducting transition temperature Tc ~ 16 K and 
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FIG. 1. Temperature dependence of thermal conductivity 
at several applied fields, plotted as k/T vs T^, for (a) H — 
0, 0.3, 1.5, 4, 6 and 8 T, and (b) H = 0, 50, 75, 100, 200 
and 300 mT, in increasing order. The solid line indicates the 
value expected from the Wiedemann- Franz law above Hc2- 



an upper critical field -ffc2(0) ~ 7 T. Its thermal conduc- 
tivity K was measured in a dilution refrigerator using a 
standard steady-state technique. A heater and two Ru02 
thermometers were used for the measurements. The later 
were calibrated in-situ for each field against a Germa- 
nium thermometer. The temperature was increased at 
fixed magnetic field, from 70 mK up, in fields ranging 
from to 8 T. The field was applied parallel to the c-axis 
of the tetragonal crystal structure ([GDI]), and perpen- 
dicular to the heat current (along [100]). Below 1.5 T, 
the sample was cooled in the field to ensure a good field 
homogeneity. Above 1.5 T, the results were independent 
of cooling procedure. The single crystal was grown by a 
melting flux method Q. The sample was a rectangular 
parallepiped of width 0.495 mm (along [001]) and thick- 
ness 0.233 mm (along [010]), with a 1.59 mm separation 
between contacts (along [100]). In zero magnetic field, 
/5(300 K) = 35 /ifJ cm and po = 1.30 fifl cm, from a 
fit to p = po + AT'^ between and 50 K. There is a 
positive magnetoresistance such that p{8 T) — 1.67 fifl 
cm at T 0. The zero-temperature coherence length 
Co = 70 A, using Hc2{0) = $o/27rCo- The zero temper- 
ature penetration depth is Aq = 760 A p5|. The mean 



free path is approximately 500 A. We measured the lower 
critical field iJci(O) to be 60 mT, using the sudden drop 
in k{H) vs at 2 K, caused by the strong scattering of 
phonons by vortices as they first enter the sample. 

The thermal conductivity k(T) of LuNi2B2C is plot- 
ted in Fig. 1, as k/T vs T^. The total conductivity 
is the sum of an electronic and a phononic contribu- 
tion: K = Ke + Kph- By plotting the data in this way, 
one can easily separate the electronic term linear in T 
from the phononic term cubic in T. As the tempera- 
ture is decreased, and in the absence of strong electron- 
phonon scattering (i.e. at _ff ~ 0), the phonon mean 
free path eventually grows to reach the size of the crys- 
tal, at which point Kph ^ T^. The low- field curves 
(in the lower panel of Fig. 1) are indeed roughly lin- 
ear (and parallel) in such a plot, with a slope in quan- 
titative agreement with the known sound velocities and 
sample dimensions, as reported earlier []l6|| . At higher 
field, roughly above 2 T and all the way into the nor- 
mal state, the conduction is essentially entirely due to 
electrons and given by a constant k/T. Note that the 
magnitude of this linear k is in perfect agreement with 
the Wiedemann-Franz law, namely k^/T = Lq/p, where 
La = (7rV3)(fci3/e)2 = 2.45 x 10"^ W fl K^^ and p = 
p{8 T). Given this well-understood behaviour of Ke{T) 
and Kph{T), it is straightforward to extract the elec- 
tronic contribution Ke{T), by simply extrapolating k/T 
to T = 0. The result of this extrapolation is plotted as 
Ke/T vs H in Figs. 2 and 3, where the field is normalized 
to unity at i?c2(0) and the conductivity, to its normal 
state value. One immediately notices the large amount of 
delocalized quasiparticles throughout the vortex state of 
LuNi2B2C. [This would seem to provide a natural expla- 
nation for the observation of de Haas- van Alphen oscilla- 
tions down to unusually low fields (i?c2/5) in YNi2B2C 
IitI , a close cousin of LuNi2B2C, with Tc = 15.5 K and 
Hc2 = 6.5 T.] In fact, the growth of quasiparticle conduc- 
tion starts right at Hd and is seen to be roughly linear in 
field. This is in dramatic contrast with the behaviour of 
quasiparticles in s-wave superconductors with a large fi- 
nite gap for all directions of electron motion. For compar- 
ison, we show in Fig. 2 the electronic conductivity of Nb 
measured at 2 K {i.e. 0.22 Tc) [^. In an isotropic s-wave 
superconductor, the only quasiparticle states present at 
T <^Tc are those associated with vortices. When vortices 
are far apart, these states are bound to the vortex core 
and are therefore localized, and unable to transport heat. 
They thus contribute to the specific heat but not to the 
thermal conductivity. As the field is increased and the 
vortices are brought closer together, tunneling between 
states on adjacent vortices will cause some delocaliza- 
tion. This conduction is expected to grow exponentially 
with the ratio of intervortex separation to vortex core 
size (~ Co) I namely as exp(— a-^/ Hc2/ H), where a is a 
constant, as is found for Nb at fields below iJc2/3 [ p^ . 

In the presence of nodes in the gap, the dominant 
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FIG. 2. Magnetic field dependence of the electronic ther- 
mal conductivity k/T at T ^ 0, normalized to its value at 
Hc2- Circles are for LuNi2B2C, squares for UPts and dia- 
monds for Nb [|l^ . Note the qualitative difference between the 
activated conductivity of s-wave superconductor Nb and the 
roughly linear growth seen in UPts, a superconductor with a 
line of nodes. The lines are a guide to the eye. 



mechanism for quasiparticle transport in the vortex state 
is totally different. Conduction results from the popula- 
tion of extended quasiparticle states in the bulk of the 
sample, outside the vortex cores. The excitation of these 
quasiparticles proceeds via the Doppler shift of their en- 
ergies as they move in the presence of the superfluid flow 
circulating around each vortex. Because near the nodes 
such states exist down to zero energy, the growth in the 
zero-energy quasiparticle density of states starts right at 
Hci, with a characteristic \/H dependence This 
leads to a dependence of the specific heat at low tem- 
perature, as observed for example in the cuprate super- 
conductor YBa2Cu307 ||2^. Note that the same mech- 
anism will operate for an anisotropic s-wave gap if the 
field is such that the Doppler shift exceeds the minimum 
gap in the quasiparticle spectrum. It is worth noting, 
however, that a similar field dependence has also been 
observed in s-wave superconductors, such as NbSe2 p^ ] 
where it has been attributed to the bound states in the 
vortex core. Specific heat studies are therefore unable to 
distinguish between a contribution coming from lo- 
calized core states and that coming from extended states 
outside the core. In contrast, thermal conductivity is 
selective, in that it only probes the contribution of dclo- 
calized excitations. 

The effect of vortices on quasiparticle transport in an 
unconventional superconductor with a line of nodes in the 
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FIG. 3. Field dependence of the electronic thermal conduc- 
tivity k/T at T — > at low fields, normalized to its value at 
Hc2- For LuNi2B2C, the growth is linear and starts at Hd, 
as emphasized by the solid line. The growth is equally rapid 
for UPta 1^ . The equivalent data for s-wave superconductor 
VaSi shows a much slower growth. 



gap function was studied in beautiful detail by Suderow 
and co-workers [ p3[ . Their measurements of k{T,H) in 
UPta yield a roughly linear increase of k^/T at T — > 
with H, shown in Figs. 2 and 3. The data is for a heat 
current in the basal plane of the hexagonal crystal struc- 
ture, which probes the equatorial line node in the gap 
function of UPts, established by transverse ultrasound 
attenuation [ p4[ . Figs. 2 and 3 reveal that quasiparticle 
conduction in the basal plane of LuNi2B2C is as good as 
in UPt^ (or even better). At low fields, the growth in the 
residual linear term kq/T = {k/T)t^o is also linear in 
H, starting at Hci: 



^ Lq H — Hci 



Pa 



c2 



(1) 



where po is the zero-field normal-state resistivity. This 
is vastly more conductive than a typical s-wave super- 
conductor. For example, electronic conduction in VsSi, 
an extreme type-II s-wave superconductor with compa- 
rable Tc (16.5 K) and (45 A), is 20 times weaker at 
H — 0.05 Hc2, as seen from data shown in Fig. 3. 

In both LuNi2B2C and UPta, the thermal conductiv- 
ity is roughly linear in H and the heat capacity follows 
approximately a \/H dependence. The latter is naturally 
understood in terms of a density of states which is linear 
in energy (coming from nodes or minima). However, a 
theory that can successfully account for the linear field 
dependence of Ke/T has not yet been formulated. 
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Nohara et al. |Tc| ] and Izawa et al. jT|] have recently 
attributed the dependence of the specific heat they 
observe in YNi2B2C to a Doppler shift of the quasipar- 
ticle spectrum as in a d-wave superconductor psf but 
applied in this case to a highly anisotropic s-wave gap, 
with a small minimum gap Amin. Interpreting the ther- 
mal conductivity data in the same way yields an estimate 
of Aniin- Indeed, because quasiparticle conduction starts 
right at Hd, the minimum gap must be smaller than the 
Doppler shift energy Eh at Hd- In a superconductor 
with a line of nodes in the gap, the average Eh is given 
by ~ Aq \/B/ Bc2 [|2^ where B is the magnetic field in- 
side the superconductor and Bc2 — Hc2- At i?ci, -B — 
in a type II superconductor thus possibly implying a true 
zero in the gap. A conservative upper bound on the min- 
imum field required to excite quasiparticles above Amin 
uses B = Hci- This gives: 



< EniHd) ^ AoV^ci/ifc2^ Ao/10 



(2) 



In other words, there is a huge gap anisotropy, with a 
minimum in the basal plane (the direction of heat cur- 
rent). 

A factor 10 in gap anisotropy is unprecedented for an 
s-wave superconductor, with a factor of 2 being the most 
that has ever been inferred in elemental superconductors 
p6t . Faced with this striking result, two questions arise: 
(1) does the gap function have s-wave symmetry (with 
deep minima) or rather another symmetry (with actual 
nodes)? (2) is the pairing indeed due to phonons? 

In relation to the first question, we stress that no siz- 
able residual linear term Kq/T is observed in iJ = (see 
p6[), a fact which would tend to argue against the pres- 
ence of nodes in the superconducting gap. However, this 
may be a question of magnitude. In the cuprates, a value 
of kq/T in excellent quantitative agreement with theory 
has been observed On the other hand, in UPta |2|] 
the observed kq /T is significantly smaller than expected 
(and in fact barely resolvable) even though there is over- 
whelming evidence for nodes. 

A possible test of the symmetry of the order parameter 
and the nature of potential nodes (imposed by symmetry 
vs accidental) is to investigate the effect of adding impuri- 
ties. While impurity scattering will reduce the anisotropy 
of an s-wave gap (either by removing nodes or by increas- 
ing Aiiiin), it will lead to more zero-energy quasiparticles 
in a gap with d-wave symmetry, for example [ p8| . Both 
Nohara et al. jl^] and Yokoya et al. have interpreted 
their data on pure and impure YNi2B2C (specific heat 
and photoemission spectroscopy, respectively) in terms 
of an anisotropic s-wave gap. 

On the question of a phonon mechanism, various au- 
thors have argued that a standard Eliashberg analysis 
of Hc2{T) and other data leads to a quantitatively satis- 
factory and consistent description of LuNi2B2C in terms 
of the measured phonon spectrum and a largely isotropic 



gap function. It remains to be seen whether such an 
analysis survives the inclusion of a very anisotropic gap. 

In conclusion, we have shown the presence of highly 
delocalized quasiparticles throughout the vortex state of 
LuNi2B2C. The quasiparticle transport grows as a func- 
tion of magnetic field in the same way as it does in 
UPta, an unconventional superconductor known to have 
a line of nodes in the gap, and not at all like s-wave 
superconductors. We conclude that the gap function of 
LuNi2B2C must have nodes in the gap, or at least deep 
minima. More work is necessary to determine precisely 
the location of these nodes/minima. Such pronounced 
gap anisotropy is unprecedented in phonon-mediated su- 
perconductivity, raising the question of whether phonons 
are indeed responsible for Cooper pairing in borocarbide 
superconductors, as has traditionally been thought. 
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